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This paper presents a second-order heavy traffic analysis of a sin- 
gle server queue that processes customers having deadlines using the 
earliest-deadline-first scheduling policy. For such systems, referred to 
as real-time queueing systems, performance is measured by the frac- 
tion of customers who meet their deadline, rather than more tradi- 
' tional performance measures, such as customer delay, queue length 

, or server utilization. To model such systems, one must keep track 

of customer lead times (the time remaining until a customer dead- 
line elapses) or equivalent information. This paper reviews the earlier 
heavy traffic analysis of such systems that provided approximations 
to the system's behavior. The main result of this paper is the de- 
^ ' velopment of a second-order analysis that gives the accuracy of the 

. approximations and the rate of convergence of the sequence of real- 

time queueing systems to its heavy traffic limit. 

O . 1. Introduction. 

— ■ 1.1. Background. Real-time queueing systems are queueing systems whose 

^ ! customers have specific timing requirements. These systems arise in voice 

and video communication systems, control systems including avionics and 
automotive, and many aspects of modern manufacturing systems. The per- 
formance measures associated with such systems relate to the ability of the 
system to meet the customers' timing requirements as a function of the 
•Oj . workload and the queue discipline. This is quite different from the more 
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common queueing system performance measures, such as customer delay, 
queue length or server utilization. In real-time queues, a queue discipline 
such as earliest deadline first (EDF) should be used rather than a more 
standard queue discipline like first-in-first-out (FIFO), which ignores the 
customer deadlines. EDF is a well studied queue discipline, especially in 
computer systems; see, for example, the monograph by Stankovic, Spuri, 
Ramamritham and Buttazzo [15]. In certain situations, EDF is optimal. For 
example, Panwar and Towsley [12] showed that, for a G/M/c queue with 
preemption, serving customers with general deadlines, EDF maximizes the 
fraction of customers that meet their deadlines within the class of work- 
conserving policies. In this paper we also treat the case of EDF with pre- 
emption. 

In real-time queues, the state space is high dimensional in that the indi- 
vidual task deadlines, the task lead times (time until the deadline elapses) or 
some equivalent information must be kept to determine whether a customer 
deadline was met when it completes its service. This high-dimensional state 
space makes exact analysis intractable; however, a heavy traffic analysis can 
be carried out. This was done by Doytchinov, Lehoczky and Shreve (DLS) 
[4] for the single node, single traffic flow case, Kruk, Lehoczky, Shreve and 
Yeung [10] for the multi-class single station case, Yeung and Lehoczky [17] 
for feed-forward networks, and by Kruk, Lehoczky, Shreve and Yeung [11] 
for multi-flow acyclic networks. Both EDF and FIFO queue disciplines were 
considered. 

In the rest of this section we summarize some of the existing results of 
the heavy traffic analysis and introduce the second-order analysis that is 
the subject of this paper. The model, assumptions and analysis are reviewed 
and made precise in Section 2 of this paper. The main results are stated in 
Section 3, and the proofs are developed in Sections 4-6. Simulations demon- 
strating the main results are presented in Section 7. 

1.2. Previous analytic results. The heavy traffic analysis of real-time 
queues begins with a sequence of queueing systems, the nth system having 
independent strictly positive interarrival times with arrival rate \( n \ and 
the customers having independent service times with mean jj^j- Assume 

the traffic intensity = = 1 — for some sequence 7^ having a 
limit 7. It follows that the scaled workload process 



where W* is a reflected Brownian motion process with drift —7 [see (2.17)]. 
To study the lead times of the customers in the queue, we introduce 



measure- valued processes Q^ n \t)(B) and W^ n \t)(B), where Q^ n \t)(B) gives 




A W(nt) 



=>W*(t) 
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the number of customers in the queue at time t having lead times in the Borel 
set B, while W^ n \t)(B) is the work at time t associated with customers in 
the queue having lead times in B. By considering B = {— oo, y], — oo < y < oo, 
one can construct the cumulative lead-time distribution of work in the queue. 
The interval (— oo,0) is of special importance because it is associated with 
work that is late. 

To characterize the limiting behavior of these measure- valued processes, it 
is convenient to define the frontier, F^ n \t), roughly the largest lead time of 
all the customers ever having received any service. Any customer with lead 
time larger than (t) has never received any service. The frontier allows us 
to divide the workload (or customers) into two parts: those customers with 
lead times not larger than F^ n \t), that is, W^(t)(-cx),F( n )(t)], and those 
with lead times greater than F^ n \t), that is, W (n) (i)(-F (n) (i),oo). DLS [4] 
prove that the scaled version ofW^ n \—oo, F^ n '(t)], namely, ^W^ n -'(nt)(— oo, 

F^(nt)], has limit zero; hence, we focus on the behavior of W^(t)(j/ V 
F^ n \t), oo). (We use the notation a V b = max{a, b} and a A b = min{a, b}.) 

Under the heavy traffic scaling implicit in the discussion so far, in order 
to obtain nontrivial limits as n — > oo , time is accelerated by the factor n and 
the workload and number of customers in queue is scaled by Because 

the unsealed workload is of order -y/n, the time each customer spends in 
queue is also of order y/n. In order to have a nontrivial limiting lead-time 
distribution, the lead times of arriving customers in the nth system must be 
of order yfn. We assume therefore that arriving customers in the nth sytem 
have lead times equal to ^Jn times random variables drawn independently 
from a cumulative distribution function G. We assume that G satisfies y* = 
min{y|G(y) = 1} < oo, so that these random variables are bounded from 
above. It is natural to assume that G(0— ) = 0, so that all lead times are 
nonnegative, but we do not need this for our analysis and, hence, do not 
assume it. 
We set 

V (1 - G(r})) drj, if y<y*, 

0, ify>y*. 
The function H maps (— oo,y*] onto [0, oo) and is strictly decreasing and 
Lipschitz continuous with Lipschitz constant 1 on (— oo, y*\. Therefore, there 
exists a continuous inverse function H^ 1 that maps [0,oo) onto (— oo,y*]. 
In [4] it is shown that, as n — > oo, the scaled frontier process 

(1.3) F^ n \t) = -^F^(nt), 



(l-G( V ))d V 



jn 

the scaled workload measure process 
(1.4) W in) (t)(B)^-^W {n) (nt)(^n~B), 
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and the scaled queue length measure process 

(1.5) & n \t)(B)±^=Q( n Xnt)(VKB) 

v n 

converge weakly to the limiting scaled frontier process 

(1.6) F*(t) = H~ 1 (W*(t)), t>0, 
the limiting scaled workload measure process 

(1.7) W*(t)(B)±[ (1-G(y))dy, 

JBn[F*(t),oa) 

and the limiting scaled queue-length measure process 

(1.8) Q*(t)(B)±(,[ (l-G(y))dy, 

JBn[F*(t),oa) 

where A = limj^oo and /i = hin^^oo ^ n "> . Since lim n _ ) . 00 p( n ) = 1, we have 
A = /z, but shall use both of these symbols to keep track of whether a term 
appears as the limit of or the limit of . This is useful to help conjec- 
ture the correct formulas if there were multiple input streams, in which case 
the limit of the sum of the arrival rates for the streams would equal /i, rather 
than the limit for any particular arrival stream. It also aids in simulation 
(see Section 7), in which we are not yet at the limit and, thus, must replace 
A and fi in certain formulas by A^™) and and these two are generally 
different. 

In summary, DLS [4] prove that, as n— > oo, 

(1.9) F*, W (n) ^W* and Q (n) Q*. 

These convergence results allow us to approximate W^ n \y,oo) for — oo < 
y < co and in terms of the scaled workload process W^ n \ We present 
those approximations and discuss their accuracy below. 

1.3. Second-order analysis. Our goal is to approximate the scaled fron- 
tier F^" 1 ' and the scaled workload measure using the scaled workload 
scalar W^ n \ and to determine the accuracy of these approximations. The 
workload process is most useful as an approximation quantity as it is the 
most easily measured. Recall lfW(t) W*(t), F^ n \t) H~ l {W*{t)), and 

oo) =► W*(t)(y, oo) = H(y V F*(t)) = H(y V H'\W*(t))). 

These suggest the following approximations for W^ n \t)(y, oo) and F^\ 

(1.10) W^(t)(y,oo)^H(y\/F^(t)), 

(1.11) W (n) (t)(y, oo) « H(y V n-\wM(t))), 

(1.12) F^ n \t)^H~ l {W {n) {t)). 
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In this paper we study the accuracy of these approximations by showing 
that the difference between the desired quantity and its proposed approxi- 
mation when scaled by n 1//4 converges to a limiting process. As n — > oo, the 
approximations in (1.10)-(1.12) become exact, and so some dilation of the 
difference of the two sides is necessary in order to obtain a nontrivial limit. 
In these approximations, t is held fixed and the process parameter y has 
been scaled by \fn [the \friB term in (1.4)— (1.5)]. This is actually a scaling 
of lead times, and when time is scaled by </n, space must be scaled by 
in order to obtain a central limit result. However, space has already been 
scaled by [the ^= term in (1.4) and (1.5)]. We must therefore multiply 

by ra 1 / 4 to partially cancel this and obtain a scaling of -^tji- Specifically, 
using [1] and [14], we prove the following results: 

(1.13) n y 4 [W( n \t)(y,oo)-H(yyF^(t))]^J*(yyF*(t)), 

(1.14) nV 4 [wW(t)(y,oo)-^(yVi7-H^ (n H0))]^^(y)Vw< S /}' 

(1.15) n^[H-\W^\t))-F^{t)] => 1 f^| t ) )) , 

where J* is a mean-zero Gaussian process with continuous paths and a 
covariance function defined in Theorem 3.3 below. In (1.13) and (1.14), both 
the left- and right-hand sides are processes in the parameter y <y* with t > 
fixed; the convergence is weak convergence in D(—oo,y*]. In (1.15), t > is 
again fixed and we have weak convergence of random variables. 



2. The model, assumptions and notation. 



2.1. The basic model. We now specify the model and its assumptions 
precisely. We have a sequence of single-station queueing systems, each serv- 
ing one class of customers. The queueing systems are indexed by superscript 
(n). The interarrival times for the customer arrival process are {v,j }^Li> 
a sequence of strictly positive, independent, identically distributed random 
variables with common mean -r 1 ^- and standard deviation a^ n \ The service 

times are {vj n ^}|2 =1 , another sequence of positive, independent, identically 
distributed random variables with common mean — ttt and standard devia- 

tion /3^ n \ We assume that each queue is empty at time zero. 
We define the customer arrival times 

k 

(2.1) s^^o, s£ n) ^X>J n) , k>l, 

i=l 

the customer arrival process 

(2.2) A^{t)=m a x{k;S { ^ <t}, t>0, 
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and the work arrival process 

(2.3) ^w=Er. 

3=1 

The work that has arrived to the queue by time t is then y( n )(^(")(t)). 

(n) 

Each customer arrives with an initial lead time Lj , the time between the 
arrival time and the deadline for completion of service for that customer. 
These initial lead times are independent and identically distributed with 
distribution given by 

(2.4) F{L<f l) < y/Hy} = G(y), 

where G is a right-continuous cumulative distribution function. We define 

(2.5) y*±mm{yeR;G(y) = l}, 

and assume that y* is finite. We assume that, for every n, the sequences 
{uj}JL 1 , {v^}'jL 1 and {Lj}^ =1 are mutually independent. 

We assume that customers are served using the earliest-deadline-first 
(EDF) queue discipline, that is, the server always serves the customer with 
the shortest lead time. Preemption occurs when a customer more urgent 
than the customer in service arrives (we assume preempt-resume) . There is 
no set up, switch-over or other type of overhead. Late customers (customers 
with negative lead times) stay in queue until served to completion. 

The netput process 

(2.6) N {n \t) = V {n) {A^(t))-t 

measures the amount of work in queue at time t, provided that the server is 
never idle up to time t. The cumulative idleness process 

(2.7) 4 - M N^(s) 

0<s<t 

gives the amount of time the server is idle, and adding this to the netput 
process, we obtain the workload process 

(2.8) W in) {t) = N^ n \t) + I {n) {t), 

which records the amount of work in the queue, taking server idleness into 
account. All the above processes are independent of the queue service disci- 
pline, provided that the server is never idle when there are customers in the 
queue. However, the queue length process Q^ n \t), which is the number of 
customers in the queue at time t, depends on the queue discipline. All these 
processes are right- continuous with left-hand limits (ROLL). 
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2.2. Heavy traffic assumptions. We assume that the following limits exist 
and are all positive: 

lim A (n) = A, lim = a, 

n— >oo n-^oo 

(2.9) 

lim a W = a, lim (5 {n) = /?. 

n-^oo n— >oo 

Define the traffic intensity p^ = . We make the heavy traffic assumption 
(2.10) lim Vn(l-P (n) )=7 

n— >oo 

for some 7 G R. This implies that X = p, although, for reasons explained 
in Section 1.3, we shall use both symbols. We use the notation p = A/ fi in 
certain formulas, even though p=l. We also impose the modified Lindeberg 
condition on the interarrival and service times: for every c > 0, 



(2.11) 



^ E [(^ > -(^ (B) )" 1 )Vr-(^)-H>«^ 

= n lim E[(^ - (M^)- 1 ) V,^,,-!,^ =°- 



This condition is satisfied, for example, if sup ngN E(ii^) 2+<5 < 00 and 

< 00 for some 5 > 0. Clearly, (2.11) implies the usual 
Lindeberg condition on the interarrival and service times: 

lim E[(*4 n) - (A (n) ) _1 ) 2 L, w a(n)WK ^] 

n— >oc LV J v ' y {\u) '-(AW) i |>c % /n} J 

(2.12) 

= lim E[(v\ n) - (/i (n) ) _1 ) 2 I (n) ]=0 

for every c > 0. It can be shown that (2.12) does not, in general, imply (2.11). 

We introduce the heavy traffic scaling for the idleness, workload and queue 
length processes 

j( n )(i) = -L/ (n) (rat)> 
(2.13) W {n) (t) = ^=W {n) (nt), 

Q(»)(t) = -LQW(„*) > 
and the centered heavy traffic scaling for the arrival processes 
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^» w ._Lg(„ r -_L), 

A^it) = -^=[A^ n \nt) - A (n W]. 

We define also 

(2. 14) N&> (t) = -4= [V (n) U (n) (rat)) - "*] • 

Vra 

Note that W^ n )(t) = iV*( n )(t) + T^(t). 

Theorem 3.1 in [13] and Theorem 7.3.2 in [16] imply that 

(2.15) (S^ n \A^)^{S*,A*), 

where A* is a Brownian motion with zero drift and variance a A 3 per unit 
time and 

(2.16) S*(\t) = --A*{t), t>0. 

A 

It is a standard result [6] that 

(2.17) (N^ n \ , W {n) ) => ( N * , I* , W* ) , 

where iV* is a Brownian motion with variance (a 2 p 2 + /3 2 )A per unit time 
and drift -7, I*(t) = - min < s < t iV*(s), and W*{t) = N*(t)+I*(t). In other 
words, W* is a reflected Brownian motion with drift —7, and J* causes the 
reflection. 

Here and elsewhere, the symbol denotes weak convergence of measures 
on the space D(T,S) of functions from a set T (which is either a closed 
interval in 1R or a closed rectangle in M 2 , both possibly unbounded) to a 
Polish space S that are right-continuous with left limits. If S = M, we shall 
write simply D(T). Throughout this paper, we shall use two topologies on 
D(T,S). In almost all places, the Skorohod J\ topology will be employed. 
This topology is convenient for dealing with weak convergence to continu- 
ous processes. The definition of the J\ topology can be found, for example, 
in [3, 5, 16] for TCI and in [1] for the case of a rectangle T in M. 2 . In 
the sequel, whenever we consider weak convergence in D(T,S), we always 
assume that this space is endowed with the J\ topology unless explicitly 
stated otherwise. In particular, (2.15) and (2.17) hold in the J\ topology. 
However, in Theorem 3.4, the (weaker) M\ topology on a half-line T will 
be used. We need to use the latter topology to establish stochastic-process 
limits with unmatched jumps in the limit process, for example, in the case 
of functions with asymptotically vanishing maximal jumps converging to an 
indicator function. See [16] for a definition of the M\ topology and more de- 
tails. We usually take T = [0, 00) and S = M. d , with appropriate dimension d 
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[e.g., d = 2 in (2.15) and d = 3 in (2.17)], unless explicitly stated otherwise. 
For stochastic processes A and B with sample paths in D(T,S), we shall 

write A = B if A and B have the same distribution on D(T, S). 

2.3. Measure-valued processes and frontiers. To study whether tasks or 
customers meet their timing requirements, one must keep track of customer 
lead times, where the lead time is the time remaining until the deadline 
elapses, that is, 



In this section we define a collection of measure- valued processes that will be 
useful in the analysis of the instantaneous lead-time profile of the customers. 
Queue length measure: 



lead time 



deadline — current time. 




Number of customers in the queue at time 
t having lead times at time t in B C R 



Workload measure: 




Work in the queue at time t associated with 
customers having lead times at time t in B C R 



Customer arrival measure: 




Number of arrivals by time t, whether 
or not still in the system at time t, 
having lead times at time t in B C R 



Workload arrival measure: 



( having lead times at time t in B C R 
The following relationships easily follow: 




(2.18) Q {n \t) = Q (n) (t)(R), W {n \t) = W (n) (t)(R) 

(2.19) A ( - n \t)=A {n) (t)(R), V ( - n \A ( - n \t)) = V {n) (t)(R), 



A^(t) 



3=1 3 3 



(2.20) 




A^(t) 



3=1 J J 
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To study the behavior of the EDF queue discipline, it is useful to keep 
track of the lead time of the customer currently in service and the largest 
lead time of all customers, whether present or departed, who have ever been 
in service. We define the frontier 

{Largest current lead time of all customers who have ever 
been in service, whether still present or not, or \/ny* — t, 
if this quantity is larger than the former one 

and the current lead time 

f Lead time of the customer in service, 



C (n) (i) 



\ or (i) if the queue is empty 



Prior to arrival of the first customer, F^ n \t) = y/ny* — t. Under the EDF 
queue discipline, there is no customer with lead time smaller than C^ n \t), 
and there has never been a customer in service whose lead time, if the 
customer were still present, would exceed F^(t). Furthermore, C^ n \t) < 
F^ n \t) for all t > 0. Both F^ and are RCLL. 

We define the scaled versions (1.3), (1.4), (1.5) and C^ n \t) = -^C^ n \nt) 
of the processes defined above under the EDF queue discpline. We define 
also 



\/n 



3=1 
^ oo 

^ £ eV^B+nt-Lf ,sf > <nt} ' 



V^(t)(B) ^^=V^ n \nt){^iB) 
\ n 



1 \ - ' (n) T 



^ oo 



E(n)| 
Jl {sj n) e>/SB+nt-£f ) > sj n) <nt}" 
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3. Main results. Denote by M the set of all finite, nonnegative measures 
on <8(R), the Borel subsets of K. Under the weak topology, M. is metriz- 
able as a complete, separable topological space. We recall Proposition 3.10 
and Theorem 3.1 of [4], which characterize the limiting distributions of the 
workload measure and the queue length measure under the EDF service 
discipline. 

Proposition 3.1 (Proposition 3.10 of [4]). We have => F* as 
oo. 

Theorem 3.2 (Theorem 3.1 of [4]). Let W* and Q* be the measure- 
valued processes defined by 

(3.1) VV*(t)(B)= [ (1-G(y))dy, Q*(t)(B) 4 f,W*(t)(B), 

JBn[F»(t),ao) 

for all Borel sets BCK. The processes and converge weakly in 

D([0,oo),M) to W* and Q* , respectively. 

By Theorem 3.2, for every t > and j/£l, 
W< n >(t)(y, oo) W*(t)(y, oo) = H(y V F*(t)) = H(y V H~ 1 (W*(t))). 
In particular, 

(3.2) VVW(t)(y,oo)&H(yVFW(t)). 

In some cases, the frontier F^ n ' (and, thus, its rescaled counterpart) may 
be difficult to evaluate. By Proposition 3.1, (1.6) and Theorem 3.2, we can 
replace F^ n \t) in (3.2) by the approximate rescaled frontier H~ l {W^ n \t)), 
getting 

(3.3) W (n) (i)(y, oo) « H(y V H'\W^(t))). 

The aim of this paper is to investigate the rate of convergence of W^(t) 
to W*(t) in Theorem 3.2. More precisely, we find the empirical processes 
corresponding to the workload measure. In what follows, we fix t > 0. Our 
main results are the following: 

Theorem 3.3. As a process in y < y* , we have the convergence 

(3.4) n l l A [W^\t){y, oo) - H(y V F^(t))} J*(y V F*(t)) 

in D(—oo,y*}, where J* is a mean-zero Gaussian process with continuous 
paths and covariance 

E[J*( yi )r(y 2 )}=a 2 Xp 2 f V f (£- yi )A(k-y 2 )dG(e)dG(k) 
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(3.5) ^_G(( yi Ay 2 ) + x)+/3 2 J 

x (l-G(( yi Vy 2 )+x))dx, 

independent of F* (t) . 

Theorem 3.4. As a process in y < y* , we have the convergence 

(3.6) n^[wW(t)(y,oo)-i/(yVi7- 1 (W(")(t)))]^J*(y)I {F , (t) <, } 

in D(—oo,y*] endowed with the Mi topology, where J* is as in Theorem 3.3. 

Let us note that Theorems 3.3 and 3.4 characterize the accuracy of the 
approximations (3.2) and (3.3), respectively. One might also ask about the 
accuracy of approximating F^ n \t) by H~ l (W^ n \t)). The answer to this 
question is given by the following proposition: 

Proposition 3.5. For a fixed t>0, we have the convergence 

(3.7) n^[H-\W^\t))-F^\t)] =► 

where J* is as in Theorem 3.3. 

Let us observe that, by (1.6), P{F*(t) < y*} = ¥{W*(t) > 0} = 1, so the 
limit in (3.7) is well defined. 

4. Customers behind the frontiers. In this section we prove that the 
work in the nth system at time nt associated with customers in this sys- 
tem having lead times smaller than the current frontier F^ n \nt) becomes 
negligible after division by n 1 / 4 . The following lemma is a refinement of the 
second part of Proposition 3.6 in [4]. 

Lemma 4.1. Under the earliest- deadline- first queue discipline, we have 

(4.1) n 1 / 4 (t) [C (n) (i) , F (n) (*)) => 0. 

PROOF. We define 

(4.2) r (n) (t) = sup{sG [0,t];(5 (n) (s)=F (n) (s)}. 

Because the system is initially empty, C^ n ' (0) = F^ n ' (0) = and the above 
supremum is not over the empty set. In the proof of Proposition 3.6 in [4], 
it is shown that 

(4.3) v^(i-r (n) (t))^0. 
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Let us observe that, by the definition (4.2) and the fact that the interarrival 
times are strictly positive, we have 

(4.4) wW(rW(i))[aW(rW(i)),F(")(rW(t))] < -L v % {nT(n){t)y 



with strict inequality only when T^ n \t) =t and C^ n \t) = F^ n \t), in which 
case the left-hand side is zero. We want to show that 



( 4 - 5 ) n l/4 V A(") (nr(«) (*)) ^ °' 

To this end, let us choose a sequence e n [ so slowly that P{A n ) — > 1 as 
n — > oo, where ^4 n = [rW (t) >t — e n j yfn\ [such a choice is possible by (4.3)] . 
By (2.15) and the differencing theorem (see, e.g., Appendix A of [4]), 

-^\A^(nt) - A^(nt - ^e n ) - \ {n) V^e n ] = A^(t) - I (n) (t - ^) => 0. 

Therefore, it is possible to find a sequence a n { so slowly that P(B n ) — > 1 
as n — ► oo, where i? n = [^4^(nt) — A^ n \nt — y/ne n ) < b n y/n], and 
b n = X^£n + o n - For any 5 > 0, we have 



lim ] 

n— >oo 



1 



CO 



(4.6) 



< lim P 

n— >oo 



< lim 1 

n— >oo 



< lim P 



n l/4' U J 4(™)(nr(™)(t)) 
1 



> (5 \ n A n n s n 



max 



X / 4 A(")(nt)-b„V^<j<A(")(nt) 3 



> 8 } n A n n 5 n 



max 



X / 4 J 4(")(nt)-fe nv ^<j<A(")(nt) J 



( n ) \ JT 



1 



max 



n 1 ^ l<i<6„VH+l 



( n ) ^ r 



lim I 

n— >oo 



n { 



(n) . 



3=1 



< 1 



lim 1 

n— >oo ' 



W / ^V4lL&nV^J+l 
„,(")\2x \b n ^H\+l 



lim P[^r ; < <5n 1/4 l 

n— >oo 



0. 



where the fourth line follows from the independence of the service times and 
the arrivals, the sixth one from the fact that {v^^JLi ls a sequence of i.i.d. 
random variables, and the last one from (2.9) and the fact that b n — > 0. But 
P(A n n B n ) -> 1, so (4.6) implies (4.5). 

An upper bound on the work with lead times in [C^ n \nt), F^ n \nt)) at 
time nt is the work arrived to the system during the time interval [nr^ (t),nt] 
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minus the work served, which is nt — nT^ n \t). Prom this and (4.4), we obtain 
the bound 



< n 1 ' 4 W {n) (t) [C (n) (t) , F (n) (t ) ) 



/ 1 O) 



+ 



1 



AW (ni) 



■;?. 



l/4"A(")(nrC")(t)) ^ -1/4 



£ ^)_ n 3/4 ( i_ r (n) (t)) 



j =J 4(")(nT( n )(t))+l 



(4.7) 



1 



n l/4 u A( n )(nr(»)(t)) 



+ n 



i/i 



y(n) ^I A (n.) (n ^^ _ y(n) Q^(")( nT (n)(t))^ 



+ ^7^[^ (B) («*) " ^ V W (*))] - n 3 / 4 (t - rW(t)) 



1 



1/4 (nr(")(i)) 



+ n 



1/4 



^W^_LlW(t) + A( n )t 



n 



+ 



r? 



i/i 



(n) 



[A( n \t) - i (n) (rW(t))] - (1 - p^)n 3/ \t - T^ n \t)). 



By (4.5), the first term on the right-hand side of (4.7) converges to zero in 
distribution. Also, by (2.10) and (4.3), 



(4.8) 



0. 



Fix e > 0. We have 

n l/4[2(™)(t)-I(™)( T W(t))] 



(4.9) 



n 



1/4 



^(n) ( t ) _ ^(n) ^ _ v r (n) 



+ n 



1/4 



2W( (t-4=) VrW(t) -i(™)(r( n )(t)) 



The second term on the right-hand side of (4.9) converges weakly to zero 
by (4.3). The first term is bounded above by 



(4.10) 



max \A^(t)-A^(s)\ 
A^ n \t)-A^ (t 



t— e / \/n<s<t 
,1/4 



max n 

0<s<e 
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One can check that ordinary and renewal functional central limit theorems 
for triangular arrays (see, e.g., [5, 7, 13] and Theorem 14.6 in [3]) imply that 



(4.11) 



!(") ( s ) 4 n 1 / 4 (AW (t) - AW (t 



B(s) 



where B is a Brownian motion (with zero drift and variance a 2 A 3 per unit 
time). Therefore, (4.10) converges weakly to maxo< s < £ |jB(s)|, which, in turn, 
converges to zero when e J. 0. We conclude that 



(4.12) 



n 



The analysis of the second term on the right-hand side of (4.7) is similar to 
that given above. For a fixed e > 0, we have 



n 



1 



n 



1 



n 



^_I(«)( t ) + \(n) t ) _ y{n) ( _^ n )(r^\t)) + A (n) r (n) (t) 

1/4 



■ n 



(4.13) 



+ n 



1/4 



_^)(^.)(( t _^) VT (-) (t) 

VtW(*))) 
VrW(t) 



+ A (n) ( (t--^= ) Vr (n) (t) 



£ 

n 

£ 

n 

e 
n 



V {n) ( ^I (n) (r (n) (t)) + A (n ¥ n) (t) 



As before, the second of the two terms on the right-hand side of (4.13) 
converges weakly to zero by (4.3). To analyze the first one, we define a 
sequence of processes (depending on the parameters t > and e > 0) by 



V t f(s)±n 



-1/4 



F (n) (y^ + A (n) n(^- 

-vw(\w n (t-A= 



By the above-mentioned functional central limit theorems for triangular ar- 
rays, 



(4.14) 



Vff(s) 



B(s), 
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where B is a Brownian motion (with drift zero and variance (3 2 per unit 
time). Moreover, 

n V4 yW|^L|»)( t ) + A («) t J 

_^(_L^,(( t __|) VT «», (t) ) 
+ a ,,(( ( _^) VtW(() ))] 

(4.15) = V${& n \t) + \V>{y/nt - {y/nt - e) + )) 

-^(^((t--L)vrW(t) 

+ A( n )(v^T( n )(t)-(v^t-e) + ) + 
+ 0(n~ 1 / 4 ). 

The right-hand side of (4.15) converges weakly to by (2.15), (4.3), (4.14) 
and the differencing theorem. This, together with (4.7), (4.5), (4.8), (4.12) 
and (4.13), gives (4.1). □ 

5. Approximation for the workload arrival measure. By Proposition 3.4 

of [4], for every y <y*, 

(5.1) sup \vW(t)(y,oo)-H(y)\^0. 

yo<y<y* 

In this section we want to find the joint limiting distribution for the the 
rescaled workload W^ n \t) and the empirical process 

(5.2) jW(y) 4 n l /\V^(t)(y, oo) - H(y)), y < y*, 
corresponding to (5.1). For y < y* , let 

M \ n \y) (B) - 

We have 



1 oo 

vw ( t )(y ,oo)=^5:^ ) i 

(5.3) 

=4 n) (y)+4 n) (y) 



{nt+0iy-£$ n) <sj n) <nt} 
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where 

rv* l 



py 1 

( 5 - 4 ) ^ = 1 ^y + V-ny-^<sr<-^ dM t ] ^ 



1 1 



In Sections 5.1 and 5.2 we study the limiting behavior of l[ n ^ and 1% , re- 
spectively. We shall see that in the limit n l ^l[ n ^ depends only on the service 

(n) 

times and lead times, not on the interarrival times, whereas obviously 
depends only on the interarrival times. Hence, the limits of n 

l/4 7 (n) 

and 

n 1 / 4 ^™'* — H), both of which exist, are independent. In Section 5.3 we use 
these obervations to characterize the limiting distribution of (J^ n \ W^ n \t)). 

5.1. Asymptotic analysis for l[ n \ For £ > y, we have 
A {n \nt) - A {n \nt - sfn~{l - y)) 

(5.6) = A (n) V^(£ -y) + Vn~ ( I (n) (t) - (t-^=(£-y) 



, n 

= A (n) v^(£ -y)+ n 1 / 4 !^ (I - y). 

From (5.4) and (5.6), by the fact that customer arrival times are independent 
of their service times and lead times, we get 



j=A( n )(nt-*Jn(l-y))+l 
* AW (nt)- AW {nt-y/n{l-y)) 



(5.7) i-Aj* E "^i B) W 

7=1 

1 / 



?1 



1/4 



r e 

^ 7=1 

For s > and y <y* , let us define a random field 

(5.8) !*>(.,,) * ^ lA '|f J («*V) - Mf + ,)). 

Then 

i L* Cn) ^MJ 
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'y 



n 



1/4 



/ E %>i/(A(»)^) +2 /}^i n) W 

A 7=1 



(5.9) 



1 LAW^Cir-yJJ, ■ 



7=1 



yW(y*-y,y-). 



Let 



^ (n)(y) -^/ £ E H 

y V 7=1 7=1 / 

(5.10) 

Then, by the second inequality in (5.7) and (5.8)-(5.10), 

(5.ii) n 1/4 /{ n) (i/) = Y(n) (y* - y,y-) + R {n) (y)- 

In the remainder of this subsection we find the limiting distribution for y( n ) 
and show that the process converges weakly to zero. 

Proposition 5.1. For every sq > andyo < y* , Y^> ^Y in D([0, so] x 
[yo,y*]), where Y is a mean-zero Gaussian random field with continuous 
paths and covariance 

E[Y( Sl , yi )Y(s 2 ,y 2 )} 

= X £ AS2 {\ G i(vl A V*) + x) + /3 2 ) (1 - G((yi V y 2 ) + x))dx. 

Proof. We will first show that the sequence {y( n )} is tight. For y <y* 
and j = 1, 2, ... , let 



Also, for < si < s 2 < s 3 < s , y < yi < y2 < V3 < y* , let 
Y^ n \(si,s i+1 } x (y k ,y k+1 ]) 

^Y^(s l+l ,y k+l ) - Y^(s i7 y k+l ) - Y^(s l+1 ,y k ) + Y^( Sl ,y k ) 



LA(")s l+lv ^J 

E 

j=LA(")s 4v ^J+l 



^74 E' Aiwf( M+1 ), 



where 



AM<»)(„, A Mf ( > + yfc) - Mf ( ' + » +1 ) 
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We define B\ = (si,s 2 ] x (2/1,2/2] and we will take B 2 to be a block "neigh- 
boring" B±, considering both the case that 

(5.12) B 2 = (s 2 ,s 3 ]x{y 1 ,y 2 \ 
is to the right of B± and also the case that 

(5.13) B 2 = ( Sl ,s 2 ] X (2/2,2/3] 

is above B\. Our goal in both cases is to obtain the bound (5.23). 

We continue with B 2 given by (5.12). The independence of the random 

variables AMj C2/Jfc 5 2/fe-t-x) for different values of j implies 

|A(")s i+ i^J 

(5.14) E(yW(B,)) 2 = T £ E(AMf ) (2/i,2/ 2 )) 2 , 

v j=LAC™)s iA /nJ+l 

(5.15) E[(y( n '(B 1 )) 2 (FW(B 2 )) 2 ] =E(y( n) ( J Bi)) 2 E(y( n) ( J B 2 )) 2 , 
for 2 = 1,2. Let us note that 

E(AMj n) (2/!,2/ 2 )) 2 



E 



(5.16) 



1 ,„ , 



: E(«< n >) 2 (GjGft) - Gjfa)) - 7-7^(^(2/2) - G s {yx)) 



<C(G j (y 2 )-G J (y 1 )), 



where 



C = supE ( ^ ) 2 = sup{(/ 3 H ) 2 + ^}< 



oc 



is a constant independent of j and n because of (2.9). Let U^ n) be a random 
variable such that 

1 



(5.17) 



Jj{n) = I 



j = l,...,[X^V^s \. 



Let X be a random variable with cumulative distribution function G, inde- 
pendent of l/( n \ For < si < s 2 < sq, yo < y± < y 2 < y* , let us define 

mp((si,s 2 ] x (2/1,2/2]) 
(5.18) =P[([/M,X -[/M) e ( S1 , S2 ] x (2/1,2/2]] 



j=Ls 1 A(") v /HJ+l 
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We have, by (5.14)-(5.18), 

E[(Y^(B 1 )) 2 (Y^(B 2 )) 2 } < c2 nA (w) ^oJ \ 2 m (n) (i?i)m (n) (j32) 

(5.19) 

<C 2 {\^s Q fmf\B l )mt\B 2 ). 

We have obtained (5.23) with = X^s^Cm^. 

Now let B 2 be given by (5.13). Using the fact that EAM^ n \y k+1 , y k ) = 
and the AMj (y k , y k +i) are independent for different values of j, we have 
El(Y<"HB 1 )f(Y<"HB 2 )f] 



n 



f L^ (n) s2v^J \ 2 / |A (7l) s 2 V^J ^ 

Vj = LA(")s lv / nJ+l / Vi=[A( n )siVnJ+l y 



[AWsav^J 

= - 2 E[(AA/f ) (y 1 ,y 2 )) 2 (AMj n) (2 /2 ,y 3 )) 2 ] 

j=LA( n )siVnJ+l 

(5.20) 

La(") S2a /hj 

+ - ^ pE(AMf \y 1 ,y 2 )) 2 E(AM k n \y 2 ,y 3 )f] 

j,fc=LA(' I )si v / nJ+l 

+ - £ E[AA/j" ) (y 1 ,y 2 )AMf ) (y 2 ,y 3 )] 

xE[AM^(y l5 y 2 )AM^(y 2 ,y 3 )]. 

Using (5.20) and proceeding as in (5.16), we can check that there exists a 
constant C\ independent of j and n such that 

E[(y( n )(Bi)) 2 (yW(B 2 )) 2 ] 

|A(")s 2v /^J 



(5.21) <-U £ (Cj(y 2 ) - Gj(yi))(Gj(ys) - Gj(y 2 )) 

n lj=LA(")s lv ^J+l 

+ 2 X! ( G i(^) " G,( yi ))(G k (y 3 ) - G k (y 2 )) . 

7,fc=|A( n )si-v/n 1+1 J 



j,fc=LAWsiVnJ+l 



Thus, 

(5.22) E[(yW(S 1 )) 2 (y^(S 2 )) 2 ] < 2( So A (n) ) 2 Cimi n) (B 1 )mi n) ( J B 2 ). 
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By (5.18), (5.19) and (5.22), for all "neighboring blocks" B x and B 2 in 
[0,s ] x[y ,y*], we have 

(5.23) E[(Y^(B 1 )) 2 {Y^(B 2 )) 2 ] < mW(Bi)mW(B 2 ), 

where m("> = X^s (C + ^/2C 1 )m ( { ) . 
It is clear that m!^ mi, and thus, 

(5.24) m (n) =>- m, 
where 

mi((ai,a 2 ] x (2/1,2/2]) ±P[(U,X - U) £ (s 1} s 2 ] x ( yi ,y 2 }}, 

(5-25) 

m = Xsq(C + \/2C\ )mi, 

with independent random variables £/ and X having distributions uniform on 
[0, so] and G, respectively. In particular, m has continuous marginals. Then, 
by (5.23), (5.24) and Theorem 3 from [1] (strictly speaking, by its extension 
described on pages 1665-1666 of that paper), the sequence {yW} is tight 111 
D([0,s o ] x [y ,y*]). (Note that vanishes along ({0} x [y ,y*])U ([0, s ] x 
{y*}), instead of the "lower boundary" ({0} x [y ,y*])U ([0,s ] x {yo}), as 
required by the assumptions of (the extension of) Theorem 3 from [1], but 
the proof of the latter result clearly goes through also in our case.) 

Now we will show that the finite-dimensional distributions of con- 
verge to the corresponding distributions of Y. Let < s\ < s 2 < sq and 
Ho < Uii U2 < y* be given. We claim that 

(5.26) l\m Q E[Y^(s 1 ,y 1 )Y^( S2 ,y 2 )]=E[Y( Sl ,y 1 )Y(s 2 ,y 2 )}. 

Let us observe that 

E[Y^( Sl ,y l )Y^(s 2 ,y 2 )]=E[Y^(s 1 ,y 1 )Y^( Sl ,y 2 )} 

+ E[Y^( Sl , yi )(Y^( S2 ,y 2 )-Y^( Sl ,y 2 ))} 

= E[Y^( Sl , yi )Y^( Sl ,y 2 )}, 

and E[Y(si,yi)Y(s 2 ,y 2 )] = E[Y(si, y±)Y(si, y 2 )] by the definition of Y, so it 
suffices to check (5.26) for s\ = s 2 = s. Assume that y± <y 2 . Then 

E[Y^(s, yi )Y^(s,y 2 )} 

= 7S E e[(<», ! ,.)-<»(^ +! , 1 )) 



x K V) -^"(a^ + »)) 
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1 [A< n >«VnJ r/ x v 

E E [(^ I 07(A Mv ^) +yi <Lf)/ v ^ } -^)( 1 -^^))j 



(5.27) 



1 V i \j/(X(")^)+y 2 <L^/^} ^(n) 



(1-G,(y 2 )) 



v - „, (n)\2. 



V" j= i L 



p^(l-G,( yi ))(l-G,( 2/2 )) 



1 



72 



— V 



(/»)' 



(I-GjCj/i)) 



(M«)) 2 
(1-^(2/2))- 



To obtain (5.26), we observe that the last term in (5.27) converges to 
A [ {^2 G ^ + x ) + ^) ( X " G ^ +*))dx = E[Y(s, yi )Y(s, y 2 )\. 

To show convergence of the finite-dimensional distributions of Y^ n \ we 
will use the Cramer-Wold device (see, e.g., [2]). Fix m, < s\ < • • • < s m < 
so, yo<yi,---,Um<y* and ti, . . .,t m € R. Then 



5> y(n) (*i,w 

m |A (n) s IV /»IJ / 



i=l 



(5.28) 



1 



72. 



1/4 



i=l j=l 

LA(") Smv ^J 

E 



5>" E X<^)-^ (b) (a5^+w 



where 



(n) a 



1 



72. 



1/1 



E *i(Mj n) (j/*)-M. 



(n)/ J 
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are independent, mean-zero random variables. By (5.28), we have 

/LA<") Smv ^J \ / m \ 

sl ^ Var E = Var E *^ (n) («. W) 



i=i 



. i=l 
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(5.29) 



so, by (5.26), 



tit fc E[yW( Si , yi )^ (n) (s fc ,y fc )], 



i,k=l 



s 2 = lim s 2 = J2 tit k E[Y{ Si , yi )Y(s k ,y k )] 



(5.30) 



i,fc=l 



Var ^2tiY(si,yi 



,i=i 



If s = 0, then E£i^ (n) (si,yi)^E^i^(si,?/i) by (5.29), (5.30) and the 
Chebyshev inequality. Assume s > 0. We shall check that the random vari- 
ables X^ satisfy the Lindeberg condition. Let = n^x!j n \ It is easy 

to see that there exist constants C±, C2 > such that |-X_- | < CiVj + C2 
for every j, n. Let e > be arbitrary. Then 



LA(") Smv ^J 

E 



1 



(X W ) 2 rf 

,1/41 3 



(5.31) < 



E 



(C?(t;j n) ) 2 + C2)dP 



SnV n ^ •/{4 n) >(VCi)(£s n n 1 /4-C 2 )} 



O(l) 



{< I, >(l/Ci)(£s n nV4_c 2 )} 



(u< n) ) 2 dP 



+ 0(1)P 



(„) SSnTl 1 ^ - C 2 



v r > 



Ci 



As 



n — ► 00, 



(5.32) 

and 
(5.33) 



( n ) eSn?! 1 / 4 - C 2 
«1 > " 



Ci 



< 



Ci 



es n n 1 /4 _ c 2 



(n) 



/r (« 



{«r>(l/Ci)(£SnnV4-C 2 )} 



(es n nV4 _ ^(n) 
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by (2.11) and (5.32). Thus, by (5.31)-(5.33), the random variables sat- 
isfy the Lindeberg condition. Therefore, by the Lindeberg central limit theo- 
rem, the sum tiY^ n \si,yi) converges weakly to iV(0, s), the distribution 
of J2iL\ tiV(si, yi). We have proved convergence of the finite-dimensional dis- 
tributions of y( n ) to the corresponding distributions of Y. 

Finally, we show that the limiting random field Y has continuous sample 
paths. By the Kolmogorov-Centsov theorem (see, e.g., [8]), it suffices to show 
that there exists a constant C such that, for any < s\ < s 2 , yi,V2 < V*, 

(5.34) E|y(si, yi ) - Y(s 2 ,y 2 )f < C\\(8 UVl ) - (s 2 ,y 2 )\\ 3 , 

where || • || denotes the Euclidean norm in K 2 . It is well known that, for every 
n, there exists a constant C n such that, for every normal random variable 
Z with mean zero, KZ 2n < C n (KZ 2 ) n . In particular, because Y is Gaussian, 
to prove (5.34), it suffices to find a constant C such that 

(5.35) E\Y( Sl , yi ) -Y(s 2 ,y 2 )\ 2 < C\\(s uyi ) - (s 2 ,y 2 )\\. 
We have 

E|y( Sl , yi )-y( S2 ,2/ 2 )| 2 

= Ey (si, yi) 2 - 2E[y(si, yi)y(s 2 , y 2 )\ + EY(s 2 ,y 2 ) 2 . 



(5.36) 
But 



|Ey(s 2 ,y 2 ) 2 -E[y( Sl , yi )y( S2 ,y 2 
= A 



(5.37) 



C G(y2 +x)+ ^) (i " G(y2 + x)) dx 

- £( K -^G((y 1 Ay 2 ) + x) + (3 2 ^ 

x {l-G{{ yi \Jy 2 ) + x))dx 

If Vi < 2/2, then the right-hand side of (5.37) equals 
A /" Sl 



+ 



-2 / 1 (G(y 2 +x)-G(y 1 +x))(l-G(y 2 + x))dx 
M Jo 

A jH [^ 2 G(y 2 +x) + /3 2 ) (1 - G(y 2 + x))dx 

^ £ ( G (y* + x ) - + x )) dx+ {^ + A/?2 ) ( S2 " Si ) 

ryz+si ryi+si\ / \ \ 

J -J )G{x)dx+^— + \(3 2 )(s 2 - Sl ) 



/J- 2 \Jy2 Jyi V/^ z 

<^ 2 \y2-yi\ + [^ 2 +\^{s 2 -s l ). 
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Similarly, if y\ > y 2 , then the right-hand side of (5.37) is equal to 

+ A jH [\G{ y2 + x) + /? 2 ) (1 - G(y 2 + x))dx 
< + /5 2 ) J*' (G( yi +x)- G(y 2 + x)) dx + \(±; + /3 2 ) (s 2 - 8l ) 
<A^+/3 2 )(|y 2 -yi| + ( S2 -si)). 

Thus, 

(5.38) \EY(s 2 ,y 2 ) 2 - E[Y(s 1 ,y 1 )Y(s 2 ,y 2 )}\<C 1 \\(s 1 ,y 1 ) - (s 2 ,y2)\\, 
where C\ = v2 A(4j + ft 2 ). By a similar (in fact, simpler) argument, we get 

(5.39) \EY( Sl , yi ) 2 -ElYisuy^Yis^y^l^dUsuyi) - (s 2 ,y 2 )\\. 
Relations (5.36), (5.38) and (5.39) give (5.35) with C = 2C X . □ 

Lemma 5.2. Let T n — o~(Uj ,j = 1, ) and let be T n -measurable 

stochastic processes such that =^ in D[0,oo). For s > and y < y* , 
let 

g (m<»>( 9 .)-m<»'(„ + ^-^)(^))). 

T/ien, /or every so > and y < y* , - Y^ =4> in D([0,s ] x [yo,y*])- 

Proof. Because -ff( n ) 0, there exists a sequence e n of positive 
numbers converging to zero so slowly that P(^4 n ) — > 1, where A n = 
[supo< s < So |-fC (n) (s)| < 

£ n \. It suffices to prove that 

(5.40) (zW-fW)I 4 ^0 

in D([0, so] x [yo,y*])- Relation (5.40) is a statement about weak conver- 
gence of stochastic processes, so the underlying probability spaces are irrel- 
evant. Thus, without loss of generality, we can assume that all the random 
variables (and, thus, all the prelimit processes) under consideration are de- 
fined on the same probability space (f2,^4,P) and, moreover, all the arrival 
times {u^}°° n=1 are independent of all the service times {wj n ' ) }^ n _ 1 and 
all the lead times {L^}j° n=1 . This is not a limiting assumption, because 
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if, for different n, the probability spaces (^ n \A^ n \f^) on which the se- 
quences {u^} ( jL 1 , {v^}JL 1 , {Lj yjLi are defined are different, we can 

take {n,A,F) = Un=i(^ {n) ^ {n) ^ {n) )- Let ? - j,n = 1, . . .)■ Note 

that, for any (deterministic) function /, 

(5.41) E[f{Y^ n \Z^)\^]=E[f(Y^\Z^)\K ( - n \-)], 

that is, we can evaluate the left-hand side of (5.41) by conditioning on a 
sample path of K^ n \-). As in the proof of Proposition 5.1, let U^ n \ n = 
1, be a sequence of random variables with distribution (5.17), let X 
be a random variable independent of this sequence and having cumulative 
distribution function G, and let (X, U^ 2 \ . . .) be independent of T . Let 

mf,"' be a random measure, depending on the sample path of K^ n \-), defined 
by 

m i 2 ) {(si,s 2 ] x (2/1, 2/2]) 

iP[([/W,I-f/W + ]fW([/W))e( Sl)S2 ] x (y u y 2 ]\K^(.)} 

where < si < s 2 < s , yo < V\ < V2 < D* and 

Let us also define random measures rri^ — 77T2 Ia„ + mi&A c , where m\ is 
the measure defined by (5.25). Observe that, for every u; £ fi, 

(5.42) m { 2 n \uj) ^mi. 
Proceeding as in the proof of Proposition 5.1, we get 

E[(zW (i?1 )) 2 (z(™)(i? 2 ))Vl^ (n) (-)] 

(5.43) <Ctmf\B 1 )M< 2 n \B 2 )l An 

<C 2 1 mi n) (B 1 )m { 2 n) (B 2 ) 

for some deterministic constant C\ and any "neighboring blocks" B±,B 2 C 
[0, s ] x [yo,y*]. As in the proof of Proposition 5.1, (5.42)-(5.43) imply that 
the random fields Z^Ia u are conditionally tight with respect to K^ n \-). 
By Proposition 5.1 and the fact that is independent of J 7 , the random 
fields (Z( n ) - Y^)lA n are also conditionally tight with respect to ifW(-). 
For any < s < sq and y$<y<y*, we have 
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E[(Z^ - Y^) 2 (s,y)l An \K W(.)] 



(5.44) 



-G[y + 



\( n h/T 



<CmP((0,s] X (y- £n ,y + en ])^0, 

where C is a constant independent of n, because, as we have seen in the 
proof of Proposition 5.1, the measures converge weakly to a con- 

tinuous measure mi. In particular, the finite-dimensional conditional dis- 
tributions of (Z^ -yW)!^ with respect to K^ n \-) converge to zero. 
Thus, the conditional distributions of the random fields (Z^ — Y^)lA n 
on jD([0, so] x [yo)2/*]) with respect to K^ n \-) converge weakly to 0, so for 
any continuous and bounded function / : D([0, sq] x [y , y*]) — > K, 

E[/((Z W _ yW)I A J|^] = E[/((ZW - Y^)I An )\K^(.)} - /(0), 

where the equality follows from (5.41). This implies (5.40), because, by the 
bounded convergence theorem, 

E/((Z (B) - yW)IaJ = E{E[/((ZW - y^)I A J| J"]} - /(0). □ 
Proposition 5.3. i?( n )=^0 in D(— oo,y*]. 

Proof. Let y < y* be given. For any 1 < j < (nt) — (nt — y/n x 
(y* - y)) and y < I < y* , A^ n \nt) - A^>(nt - y/n(l - y)) > j if and only 
if nt — \fn{l — y) < ^ A ( n )f nt \j +1 which, in turn, is equivalent to £ > y + 

i(»*-sS„„«_,x,)-Thu S , 



r 

-hi 



A( n )(nt)-j+\> 
* A^>{nt)-A^>{nt-^/n{i-y)) 

X; dM { p\i) 

3=1 



A( n ){nt)-A( n ){nt-yfti{y*-y)) . 
(5.45) = £ {Mf{y* 

j=i V 



_M 7 (n) ( y + -L(nt-S W 



A( n )(nt)~j+1- 



For s > and y <y* , let 
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where, for u > 0, 

4 ^[^» )(nt)+1 _ AWvfeJV0 " (nt - v^)] 

— + 

Observe that, by (2.15), (2.16) and the differencing theorem, with probabil- 
ity approaching 1, 

H[n) ^) = sg )t+o(1) + ^ {n \t) + o(-±=) =► o. 

By Lemma 5.2, for every so > and yo < y*, 

(5.46) xW - y< n ) =► 

in D([0,sq] x [yojy*]) and, thus, by Proposition 5.1, 

(5.47) => Y 

in D([0,s ] x [y ,y*],R). By (5.10), (5.45), the definition of X^ n \ (5.6) and 
(5.9), 

RW{y) = I'"' (y* - y + -^_l^(y* - y),y) - Y^(y* - y,y-) 

= (xw (y* - y + 3^r7i^ (n) (if - y) , v) - * (n) (y* - v, v)) 

+ (X^(y*-y,y)-Y^\y*-y,y)) 

+ (Y^(y* - y ,y)- Y W(y* -y,y-)), 

which converges to in D[yo,y*] for every yo < V* by (4.11), Proposition 
5.1, (5.46), (5.47), continuity of the sample paths of Y and the differencing 
theorem. □ 

5.2. Asymptotic analysis for . The analysis of the limiting behavior 
of n 1 / 4 /^ (y) is easier. We have the following: 

LEMMA 5.4. We have n 1 / 4 ^ — H) =>• Z in D( — oo,y*], where Z is a 
mean-zero Gaussian process with continuous paths and covariance 

E[Z( yi )Z(y 2 )] = Xa 2 p 2 f f (£ - yi ) A (k - y 2 ) dG(l) dG(k). 

Jyi Jyi 
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Proof. For y < y*, let Z(y) = ± jf B{£ - y) dG(£), where B is a zero- 
drift Brownian motion with variance a A 3 per unit time. It is easy to see 
that Z is a mean-zero Gaussian process with continuous sample paths. Fur- 
thermore, 



E[Z( yi )Z(y 2 )]=E 



f y B(£- yi )dG(£)-- f V B(k-y 2 )dG(k) 
E[B(£ - yi)B(k - y 2 )] dG{£) dG{k) 



a 2 i 

M J yi J V2 



rv cv 
„2\„2 1 



a 2 Xp 2 / (£- yi )A(k-y 2 )dG(£)dG(k). 
Jyi 'yi 



Finally, by (5.5), 



I^\y) = ^—( V (A^(nt)-A( n \nt-V^(£-y)))dG(£) 
(5.48) =^)J y V ^ n \£-y) + A^(t)-A^(t- e -^J^dG(£) 

= pW J* (£-y)dG(£) + j* - ( ± IJI AW(£-y)dG(£), 
and thus, by (2.9), (2.10), (4.11) and the fact that 



(5.49) / (e-y)dG(£) = H(y), 

we have [with A defined by (4.11)] 
n 1 /\l { 2 n \y)-H{y)) 

(5.50) 

f ^ {n) (< - *) dG W + (^74) =► □ 

5.3. Asymptotic analysis for (J^ n \W^ n \t)). 

Corollary 5.5. We have J (n) =^ J* in D(-oo,y*], where J* is a mean- 
zero Gaussian process with continuous paths and covariance (3.5). 

Proof. By (5.7), (5.9), (5.10), (5.50), Propositions 5.1, 5.3 and the 
independence of the arrivals, the service times and the lead times, 

(n 1 / i l[ n \y),n^(li n \y)-H(y))) 
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/ -, . y * A(") v ^(^-y)+n 1 /4A(")(£-y)-l 



-LiW(^- y )dG(^)J+ (l) 

aW v^(<-y)+n 1/4 -A (n) (t-v) 

(5.51) ^(^y; e ^ b) w, 

" -^-A^(l-y)dG(t)]+o(l) 

Y^(y*-y,y-) + R( n \y),J* -Li™ (£ - y) dG(t)j + o(l) 

^(Y(y*-y,y),Z(y)), 

where Y and Z are as in Proposition 5.1 and Lemma 5.4, independent of 
each other. Thus, by (5.2), (5.3) and (5.51), J*, as claimed. □ 

Proposition 5.6. We have (jW,wW(t)) => {J*,W*{t)) inD(-oo,y*} x 
K, where J* is as in Corollary 5.5 and VF* is a reflected Brownian motion 
with variance (a 2 p 2 + /3 2 )A per unit izme and drift —j, independent of J* . 

Proof. Fix an arbitrary yo < y* . By (2.17), Corollary 5.5 and its proof, 
it suffices to show that W^ n \t) is asymptotically independent of the pair 

of processes (n l / i I i ( l \y),n l / i (I l ^ l \y) — H(y))), yo < y < y* . We assume 
throughout the proof that n is sufficiently large so that nt — \/n(y* — yo) > 0. 
We note at the outset that, by (2.17) and the differencing theorem, 

(5.52) (t " - W - yo)) - W^(t) 0. 



Denntp T (n) A <? (n) and T (n) A <? (n) sn that 

Denote 1 1 - ^(n)^^^.^)) and J 2 - a A<»)(nt-VfI(ir-ito))+i' so timt 
T 1 (n) < nt - vW - yo) < T 2 (n) . Define 0< n ) = T 2 (n) - (nt - ^/E(y* - yo)). The 
process A^ n \s), < s < nt — y/n(y* — yo), and the process 

A in \y) = A^(T^ + ^(y - yo)) ~ A^{T^) 

= A&) (nt - Vn~(y* ~y) + {n) ) ~ A {n) (nt - ^(y* - y ) + # (n) ) , 

Ho < V < y*, are independent. 
We show now that 

(5.53) ^ (n) ^°" 
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Using the process that has the continuous limit S* in (2.15), we may 
write 



max |5 (n) (r)-S (n) (r-)| = max 

0<r<2At Jn l<j<\_2\nt\ 



{n) 1 

Ua 



3 \H 



. 1 (n) 1 

> max —=u) 



i<j<\2\nt\ yjn 3 s/nX&y 

Because the limit of S 1 ^ is continuous, the term maxK^x^AntJ ~J^ u j con ~ 

verges to zero in probability. Since =^ A* , we can choose a sequence of 
sets {5 n }£° =1 with F(B n ) -> 1 such that ,4( n )(nt) < 2A?it - 1 on 5 n . We set 

S W = 1 („( _ r 1 (») )> S W = 1 (nt _ T W), 



so 

_ Cri n ) rp{n)^ _ (n) (n) . 1 fl ( n ) 
^ U A(»)(nt-Vn( W '-i«))+l - - J l )~ S 1 ~ S 2 ■ 

In particular, (s^ — s^^Bn converges to zero in probability, and hence, so 
does Sj 1 — s 2 . We have s 2 = 4 — (sj — s 2 ) >y* — yo — (s^ — s 2 )' so 
IP(Cn) — ► I, where C n = [s 2 > 0]. On the set C n , the differencing theorem 
now implies that zero is the limit in probability of 

i^(4 n) ) - l (n) (4 n) ) = -^(A^(nt - 4 n) y^) - A^(nt -s^V^) 



n 



A(™)^S n) -4 n) )) 



= ^74(l-AW^(4 n) -4 n) ))- 

We conclude that n 1//4 (s^ — s^) =^0. This implies (5.53). 
Recall that 

* A( n )(nt) 

™ 1/4 4%) = ^74 f E ^i n) W. ^ < y < i/*. 

y j=A(")(nt-v^(^-y))+l 

We define the related process l[ n \y) by 



y <j=A( n )(7it— \/n(i-i/)+0( n ))+l 



71 



1 r(y+*<»>A/S)A,,* A(»)(nt+«W) 

^ j=A( n )(nt- % Ai(^-y)+f ( ™ ) )+i 



1/4 
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(5.54) 
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+ 



n 



+ n 



1/1 



AW (nt+dW ) 



/4j(n) 



(»+«W/VS)AW_. =A(B)(ftt)+1 

6>(n) 



dMj n) (£) 



Ay- 



Recall further that, by (5.48) and (5.49), 

,1/4 



nV4 ( /(%) _ F( y) ) = !L_ £ [>>(*) - (t - ^) 



+ n 1/4 (p (ri) -l)#(y), y <y<y*- 



We define the related process (y) 



n 1 /*(j<%)-.ff( y )) 



1/4 , y * 



A™ [ t + 



6>W 



+ 



dG(£) 



n 



1/4 /•( J /+e(")/yH)Ay* r 



n( n ) i 
^ I Jy 



( t + 



0{n) 



n 



(5.55) 



P f t 



re 



V4//H 



y + 



n 



Ay 
6>M 



+ 

n 



11 



dG(l) 



H y + 



+ 



n 



i/i 



(n) 



0(n) 



G(y*)-G(y + 



Q{n) 



Note in (5.55) that G{y) = G{y*) = 1 and H(y) = H(y*) = for y > y* , so 
evaluating G and H at y + gives the same result as evaluating these 

functions at (y + ^j=) A y*. The pair of processes (l[ n \y) , I^iy); yo < y < 

y*) is independent of the random variable W^ n \t — -^{y* — yo)), a s we 

now explain. In l[ n \y), the sum j = A^ n \nt — ^Jriit — y) + 9^>) + 1 to 
A^ n \nt + 6>W) has 

^(») (n£ + W ) - (nt - - y) + (n >) = A H (y* ) - A (n) (y*-(l-y)) 
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terms, and this is independent oiW^ n '(t— -^(y* — yo))- The Mj processes 

(n) 

appearing in this sum involve service times v - , but these particular indices 
j do not appear in the definition of W^ n \t — -^=(y* — yo))- The integrand 
appearing in l\ (y), 



t 



+ 



n 



A^[t-^ + 



n n 



[A^ (nt + 6 {n) )-A^ (nt - - y) + 6> (n) ) - A (n) v^(Z - y )] 
-UA (n V) - A (n \y* -(I- y)) - - y)], 



is also independent of W^ n \t — -^=(y* — yo))- 
It remains to show 



(5.56) nV\iP(y)-l[ n \y))^0, n^(jW (y ) _ . 

This will imply that the limit of the pair (n 1 ^ 4 I^ n \y),n 1 ^ 4 (I^ l \y) — H(y))) 
is independent of W*(t), the limit of W^ n \t - -j^(y* - y )). Since J (n) (y) = 

n 1/4 ^ n) (!/) + n 1/4 (4 B) (!/) ~ H (v)) J *0/)' we wil1 have the desired result. 
From (5.54), we have 



1/4 /f(n). 



r(«)/ 



n 1/4 (i? B) (i/)-/l B) (y)) 



n 



1/4 



(5.57) 



A( n )(nt+0( n )) 

^ cZM,f n) ^) 

j=A(") {nt—Jn(l-y)+e<- n ) )+l 



At") (rat+0(«)) 



n 



i/i 



j=A(")(nt)+l 



+ n 



1/4 



61 W 



Ay* 



4 n \y) 



For every y G [yo,y*], the absolute value of each of the first two terms on the 
right-hand side of (5.57) is bounded above by 



(5.58) 



1 



)l 



1/4 



E 



j=A( n )(nt)+l 



(n) , 1 
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The ordinary and renewal functional central limit theorems for triangular 
arrays (see, e.g., [5, 7, 13] and Theorem 14.6 in [3]) imply that 



(5.59) 



C< n > (s) = -Lr [A {n) (nt + V^s)- (nt) - A (n) ^s] 
n 1 / 4 



n 



1/4 



a» I t + 



A( n \t) 



C*(s), 



where C* is a Brownian motion with zero drift and variance a A 3 per unit 
time. In particular, C <Jl \sn~ 1 l 4: ) in D[0,oo) and, thus, 



(5.60) 



4tj [A (n) (nt + sn 1/4 ) - A™ (nt)} Xs 
n 1 / 4 



in D[0, oo). Let e > be given. The convergence (5.53) implies that ¥(D n ) — ► 
1, where £> n = [6»W < en 1 / 4 ]. On L> n , the expression (5.58) is dominated by 



A^int+en 1 / 4 ) 



n 



1/4 



E 



«< n) + — 
3 /i( n ) 



j=yl( n )(nt)+l 

+ en 1//4 ) — A*" - ) (nt) 



n 



1/4 



AW (ni + en 1 / 4 ) - A( n )(nt) 



^"'(nt+en 1 / 4 ) 

E 

j=A( n )(nt)+l 



,(») J_ 



which converges weakly to Xe ■ f = 2e by (5.60) and the law of large numbers 
for triangular arrays, together with the independence of the arrival times and 
the service times. Thus, the expression (5.58) converges to zero in probabil- 
ity. To see that the third term on the right-hand side of (5.57) converges to 

zero, we use (5.53), the fact that n 1 / 4 /^ (y) =>■ Y(y* — y,y) by (5.11) and 
Propositions 5.1, 5.3, the joint continuity of Y and the differencing theorem. 
This concludes the proof of the first convergence claimed in (5.56). 
For the second convergence claimed in (5.56), we use (5.55) to write 



(5.61) 



n 1 '\it\y)-i { 2 ) {y)) 

n l/4 r (y+e( n ) A/H)Ay* 



/' 



(n) 



A» [ t + 



Q{n) 



n 



A» t 



n 



n 



dG{£) 



+ n 



1/4 



r(n) 



y + 



0(n) 



Ay* 
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H y 



0(n) 



(4 n \y) - H(y)) 



n 



+ 



n 



1/4 



(») 



A™ ( t 



+ 



n 



G{y*)-G[y + 



The absolute value of the first term on the right-hand side of (5.61) is 
bounded uniformly in y Q < y < y* by sup o < s < 0( „ )/v ^ |C (n) (s)| , which 
converges to zero by (5.53) and (5.59). Convergence to zero of the second 
term on the right-hand side of (5.61) follows from Lemma 5.4, (5.53) and 
the differencing theorem. The absolute value of the third term in (5.61) is 
bounded uniformly in yo <y <y* by n l ^\p <yn ^ > — l\H(yo), which converges 
to zero by (2.10). Finally, the last term on the right-hand side of (5.61) 
converges weakly to zero by (5.53), (5.59) and the differencing theorem. □ 

Corollary 5.7. We have (jW, W^ n \t),F^ n \t)) (J*,W*(t),F*(t)) 
in the space D(—oo,y*] x M 2 . 

Proof. By the definition of the frontier, customers with lead times at 
time nt exceeding F^ n '{nt) have not received any service by that time. Thus, 
by Corollary 3.8 of [4] and the display above (3.32) in [4], 

W<ri(t) = V (n) (t)(F (n) (t), oo) + W (n) (t)[C^\t),F^(t)} 

(5.62) =V^(t)(F^(t),oo) + o(l) 

= H(F^ n \t)) + o(l). 

The function H -1 is continuous, hence, uniformly continuous on each bounded 
interval [0, c] , c > 0, and the Lipschitz constant of H~ l on [c, oo) approaches 
1 as c— ► oo. In particular, H^ 1 is uniformly continuous on [0, oo), so (5.62) 
implies that iJ^O^O^i)) - F^it) 0. Thus, by Proposition 5.6, 

(j( n \W^(t),F^(t)) = (j( n ),W( n )(t),if- 1 (t?( n )(i))) + o(l) 

=>(J*,W*(t),F*(t)). □ 

6. Proofs of the main results. In this section we prove Theorems 3.3, 
3.4 and Proposition 3.5. We need the following refinement of Lemma 4.1 to 
include the right endpoint F^ n \t). 



Lemma 6.1. We have nV4yy(»)(i)[c( n )(*),.F( n )(*)] 0. 
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Proof. By Proposition 3.1 and Lemma 4.1, it suffices to show that, for 
any fixed y <y*, 

(6-1) ^W^WII^^^^O. 
Take a sequence e n [ such that e n = c^tiT 1 / 4 ). Then 

< n V^)(^(»)(t)}i {?W( ^ } 

< n 1 / 4 [VW (t) (FW (t) - e n , oo) - V<»> (t) (t) , oo)]!- 

(6.2) 

< [J»(FW(f) - ^ - J^(F^(t))]I {9(n)myo} 

+ n l l\H(F^ n \t) - E n ) - H(F^(t))) 

<w(.?W e n )+n 1/4 £„, 

where, for x G D(— oo,y*] and (5 > 0, 

u;(x,<5)= sup \x(s2) — x(si)\ 
yo-ei<s 1 <s 2 <y* 
S2—si<6 

[and we have used Lipschitz continuity of H in the last line of (6.2)]. The 
right-hand side of the last inequality in (6.2) converges weakly to zero by 
Corollary 5.5 and the choice of e n . This shows (6.1). □ 

Proof of Theorem 3.3. For every y <y*, we have 
n 1 / 4 [W^(t)(y,oo)-H(yVF^(t))} 

= n 1 / 4 [W ( - n \t)([C^(t),F^(t))n(y,oo)) 

(6.3) 

+ V (n) (t) (y V F&> (t) , oo) - H(y V F (n) )] 

= j(n) (yVj pW(t)) + (l), 

where the first equality follows from the definition of the frontier as in the 
proof of Corollary 5.7 and the second one from Lemma 6.1. The mapping 
<& :D(— oo,y*] x (— oo,y*] — > D(-oo,y*] defined by $(x, y)(-) = x(y V •) is 
continuous on C(— oo,y*] x (— oo,y*]. Thus, by (6.3) and Corollary 5.7, 

n V A ffi n \t)(y,oo)-H(yVF<ri)] = <^>(J^ ,F^(t))(y) + o(l) 

=>*(J*>**(*))G/) 
= J*(yVF*(t)) 

in £>(-oo,y*]. □ 
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PROOF of Theorem 3.4. Let ^ : [0, oo) 2 xi-> D(— oo, y*] be the map- 
ping defined by *(o, a, c)(y) = d { . y<H -i (a)} , *(a,6,c)(y) = 7 ^ B {aAH(y)-bA 
H(y)) for a 7^ 6. It is easy to check that ^ is a continuous map of [0, oo) 2 x R 
into the space D(— 00, y*] endowed with the M\ topology. As in (6.3), we 
have 

n 1 l i [W^ n \t)-H{F^ n \t)))=n 1 / i [W {n \t)[C^ n \t),F^ n \t)) 

+ V (n) (t)(F (n) (i),oo) - H{F^ n \t))} 
= J (n) (F (n) (t)) + o(l), 

so, by Corollary 5.7, we get 

{w( n Xt),H{F( n \t))y/*[wW(t) - H(F^(t))}) 

(6 ' 4) ^(W*(t),W*(t),J*(F*(t))), 

because the mapping &:D(— 00, y*] x (— 00, y*] — ► R defined by $>(x,y) = 
x(y) is continuous on C(— 00, y*] x (— 00, y*]. By (1.6) and (6.4), 

n 1/4 [F(y V H^iW^Xt))) - H(y V F (n) (t))] 

= n 1 / 4 [(H(y) AWW(t)) - (H(y) A H(F^(t)))} 

(6.5) = ^(fW^.ff^W^)),^ 4 ^^^ - H(F^ n \t))])(y) 
^V(W*(t),W*(t),J*(F*(t)))(y) 

= J*(F*(t))I {y<F * m , 

where weak convergence holds in the M\ topology. The J\ topology on 
D(— 00, y*] is stronger than the Mi topology, so, by Theorem 3.3, (3.4) 
holds in the M\ topology. It is clear from the above argument and the proof 
of Theorem 3.3 that the convergence (3.4) and (6.5) is, in fact, joint [because 
both (3.4) and (6.5) follow from Corollary 5.7 by the continuous mapping 
theorem]. Thus, 

n 1 / 4 [yy(«)(t)(y, 00) - H{y V H^iW^it)))) 
= n 1/4 [W (n) (t)(y, 00) - H(y V F (n) (*))] 

- n 1/4 [# (y V H-^W^it))) - H(y V F (n) (t))] 
^J*(yVF*(t))-J*(F*(t))I {y<F * m 

= J*(y) I {F*(t)<y} 
in the Mi topology. □ 

Proof of Proposition 3.5. Recall from (6.4) that 

(6.6) n 1/4 [W (n) (t) - H(F^ n \t))} J*(F*(t)). 
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We derive (3.7) from (6.6) using the delta method. Let us observe that H is a 
convex, decreasing function on R. In particular, H" 1 is convex on [0, oo) and 
both H and H~ l have one-sided derivatives at each point of their domains. 
Clearly, H'(y+) = G(y) - 1, y € R. Also, 

(6.7) (H- 1 Y(H(y)-)H'(y+) = 1, y<y*. 
Indeed, for every e > and y <y* , we have 

1 =}-{H-\H(y + e))-H-\H{ y ))) 

H-\H{y + e))-H-\H{y)) H(y + e)-H(y) 
H(y + e)-H(y) e 

- (H- i y(H(y)-)H'(y+) 

as e I 0, so (6.7) holds. For every n, there is an element (depending 
on the elementary event cu) belonging to the sub differential of H~ l at some 
intermediate point z £ [W^{t)AH(F^ n \t)),W^ n \t)\/H(F^ n \t))} such that 

(6.8) H- l {W^ l \t)) - F (n \t) = D^[W {n) {t) - H{F^ n \t))]. 
By (6.7), for any z > 0, 

(6.9) (H^y(z-) 



H'(H- 1 (z)+) G(H~ 1 (z)) — 1' 
Applying (6.9) to a sequence z n [ z, we get 

(6.10) (H~ 1 )'(z+" 



G{H-\z)-)-l 
Using (6.9), (6.10) and convexity of H , we get 

1 

Gfff-^WW A tf(F(") (£)))) - 1 

(6.11) 

< flW < 



Glff-H^Ht) Vi?(F(»)(t)))-) - 1 



Multiplying both sides of (6.8) by ra 1 / 4 , using (6.4), (6.11), together with 
continuity of the distribution of F*(t) and the fact that G has at most 
countably many jumps, we obtain 

(6.12) n^[H-\W^(t)) - F»)(t)] =* cP^frT 

Finally, let us observe that J* is a mean-zero Gaussian process, indepen- 
dent of F*(t), so (J*,F*(t)) and (-J*,F*(t)) have the same distribution. 
Therefore, (6.12) implies (3.7). □ 
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7. First-in-first-out simulations. To illustrate our results, we consider 
the special case of constant initial lead times, that is, 

(7.1) G(y) = I{ y >y*} for some y* £ K. 

In this case, all customers in the nth system arrive with initial lead time 
\priy* and EDF reduces to the well-known first-in-first-out (FIFO) service 
discipline. By (7.1), we have H(y) = (y* — y) + and H^ 1 {w) = y* — w for 
all w > 0. Thus, by (1.6), F*(t) = y* - W*(t). Evaluating the right-hand 
side of (3.5), we get E[J*( yi ) J*{y 2 )] = \{a 2 p 2 + (3 2 )((y* - y x ) A (y* - y 2 )), 
so J*{y) = B(y* — y), y <y* , where B is a Brownian motion with zero drift 
and variance \(a 2 p 2 + /3 2 ) per unit time, independent of W*(t). By Theorem 
3.3, 

(7.2) n^iW^my, oo) - (y* - (y V F (n) (t)))] B((y* - y) A W*(t)) 
in D(—oo,y*] [by definition, F^ n \t) < y*]. Theorem 3.4 yields 

(7.3) n l / 4 [W^(t)(y, oo) - ((y* -y) A ^W(t))] 5(y* - I/) V(*)>iT-v} 

in -D(— oo,y*] endowed with the Mi topology. Finally, because W*(t) > 
almost surely, we have G(F*(t)) = almost surely and by Proposition 3.5, 

(7.4) n^iiy* - W {n) (t)) - F^{t)} B(W*(t)) 
in D( — oo, y*\. 

This case (with y* = for convenience) was already considered in [9] . In 
particular, (7.2) and (7.3) are contained in Corollary 4.5 and Theorem 4.1 
from [9], respectively. Note that the convergence in Theorem 4.1 from [9] 
takes place in D(— oo,0] endowed with the M\ topology, although this was 
not written explicitly in the statement of that theorem. 

We conduct two types of simulations. For the first simulation, we observe 
that an important use of the theory in [4] is to predict the amount of late- 
ness that will be incurred by the queueing system, and the result of this 
paper concerns the accuracy of that prediction. In particular, [4] provides 
the limiting lateness result (see Theorem 3.2) 

(7.5) VV*(i)(-oo,0] = W*(t) -W*(t)(0,oo) = (W*(t) -y*) + . 
In particular, 

(7.6) (W {n \t) - y*) + - >V (n) (t)(-oo, 0] 0. 

The first term on the left-hand side of (7.6) is the lateness predicted by 
the theory as a function of the scaled workload, and the second term is the 
actual lateness. Setting y = in (7.3), we obtain 

(7.7) n l /%m n \t) - y*) + - W^(t)(-oc, 0]] => B{y*)I {w * {t) > r} . 
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With the change of variable s = nt, we conclude from (7.7) that 

(W {n) (s) - v^y*) + - W (n) (s)(-oo,0] 

= Vn(W( n \t) -y*) + - VnW (n) (t)(-oo,0] 
(7.8) &n l / A B{y*)l {w * {t) > r} 

= B (Vny*) I { y ^W*(t)>^iy*} 

where W(s) = y/nW*(^) is a reflected Brownian motion with drift — 

and variance \(a 2 p 2 + f3 2 ) per unit time. 

We assume in this section that 7 is strictly positive. Under this assump- 
tion, the stationary distribution of W is exponential with mean 3, where 
6 = ^x(a^p i +/3 i ) " When conducting simulations on a single system, as op- 
posed to a sequence of systems indexed by n, we know only the prelimit 
quantities = 1 — p^ n \ X^ n \ and (3^ for a single value of n, and we 
do not know the value of n. We thus define 

0(n)A 2(1 



AW[(aW) 2 (pW) 2 + (/?W) 5 



a quantity we use as a surrogate for 9. In particular, an approximate density 
for W(s) is 

(7.9) f(y) = 9^exp(-e^y), y > 0, 

and an approximate distribution for (W^(s) — \/ny*) + — W^ n \s)(— 00, 0] is 

a mixture of a normal distribution with mean zero and variance 2 ^ 1 ~ P e( l) j ^™ v 

with total mass exp(—9^y/ny*) and a point mass of size 1 — exp(— 6^ n 'y/ny*) 
at 0. We will shortly present simulations to assess the accuracy of this ap- 
proximate distribution. 

While (7.2) and (7.3) are assertions about convergence of processes, (7.4) 
is an assertion about convergence of random variables. It describes the ac- 
curacy of the prediction of the frontier as the function H~ l applied to the 
workload. From (7.4) we have 



f^ y * _ w {n \s) - F (n) (s) = y^y* - yfHW [n \t) - y^F (n) (t) 
(7.10) &n l/i B(W*{t)) 

= B(y/nW*(t)) = B(W{s)). 
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Because B(y) is normal with mean-zero and variance approximately equal 

to 2{1 ~ g C ))V and B (y) is independent of W*(t), the density of B{W{s)) is 
approximately 



Q(n) ,oc l / Q(n) x 



2 



(7 - n) 9{x) = v 4.(1 - P H) loo ^ exp Ud - pwjfM d y- 

Using the Laplace transform formula 

e -(am-PV dy= l e -V*P ) a >0,p>0, 

o y/y VP 

we can simplify (7.11) to obtain Laplace's density 

(7.12) g(x)= =exp 1 x£R. 

2A/l-p( n ) v Jl-p( n ) / 



To test the predictive value of (7.8) and (7.10) as approximations to the 
empirical distributions of (W n (s) — y/ny*) + — W^ n \s)(— oo, 0] and y/ny* — 
we simulated a single server queueing system with Poisson 
arrivals (A^ n ) = 0.96) and with three different service distributions: Exponen- 
tial^), Gamma(2, 0.5) and Uniform[0.5, 1.5]. Each of these distributions has 
mean 1 but different variances, so in all cases, =0.96. All customers had 
a constant initial lead time of ^/ny* = 30. The queueing system began in 
an empty state and was simulated for T = 4,000 time units to ensure that 
the equilibrium assumption underlying the use of (7.9) as an approximate 
density for W(s) was valid. The simulation was independently repeated a 
total of 4,000 times to determine the empirical distribution. 

Figures 1, 2 and 3 present normal Q-Q plots of the values for (W^ n '(s) — 
y/ny*) + - W (n) (s)(-oo,0] restricted to the situation in which W^ n \s) > 30. 
In this case, these values should be approximated by a normal distribu- 
tion and the normal plot should be close to linear. These three figures 
provide strong confirmation of the accuracy of the continuous part of the 
proposed limiting distribution. Table 1 compares the theoretical probability 
that W^ n \s) > 30 with the empirical probability derived from the simula- 
tion. The tabled probabilities are the total mass associated with the con- 
tinuous part of the distribution. Again, the table shows the theory to be 
remarkably accurate. 

The same simulations were used to assess the second limiting result. 
In particular, Figures 4, 5 and 6 address the accuracy of approximating 
V^j/* - W^ n \s) - F(")(s) by the Laplace distribution defined by (7.12). 
The figures differ only in the particular choice of service distribution. The 
figures on the left present the sorted values of \Jny* — W^ n \s) — F( n \s) 
from 4,000 independent simulations of an M/G/l queue stopped at time 
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Fig. 2. 



normal plot of continuous part of (7-8): gamma service. 





Fig. 4. Sorted values of (7-10) and Q-Q plot versus Laplace distribution (7.12): expo- 
nential service. 
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Fig. 6. Sorted values of (7.12) and Q-Q plot versus Laplace disribution (7.12): uniform 
service. 



T = 4,000. To evaluate the appropriateness of the proposed Laplace dis- 
tribution approximation, the plots on the right present a Q-Q plot with 
respect to the Laplace distribution. The more linear the plot, the more ap- 
propriate is the Laplace distribution. The curves in each of these figures are 
highly linear; hence, they offer strong evidence of the appropriateness of this 
approximation . 
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Table 1 

Mass of the continuous part of (7-8) 



Service distribution 


Variance 


Theory: exp(-3O0( n ') 


Simulation 


Exponential 


1.000 


0.2865 


0.2840 


Gamma(2, 0.5) 


0.500 


0.1889 


0.1880 


Uniform 


0.083 


0.0995 


0.0975 
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